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Abstract 

We study a system of ./V non-intersecting Brownian motions on a line seg- 
ment [0, L] with periodic, absorbing and reflecting boundary conditions. We 
show that the normalized reunion probabilities of these Brownian motions in 
the three models can be mapped to the partition function of two-dimensional 
continuum Yang-Mills theory on a sphere respectively with gauge groups 
U(iV), Sp(2iV) and SO(2iV). Consequently, we show that in each of these 
Brownian motion models, as one varies the system size L, a third order 
phase transition occurs at a critical value L = L C (N) ~ y/N in the large 
TV limit. Close to the critical point, the reunion probability, properly cen- 
tered and scaled, is identical to the Tracy- Widom distribution describing the 
probability distribution of the largest eigenvalue of a random matrix. For 
the periodic case we obtain the Tracy- Widom distribution corresponding to 
the GUE random matrices, while for the absorbing and reflecting cases we 
get the Tracy- Widom distribution corresponding to GOE random matrices. 
In the absorbing case, the reunion probability is also identified as the maxi- 
mal height of N non-intersecting Brownian excursions ( "watermelons" with 
a wall) whose distribution in the asymptotic scaling limit is then described 
by GOE Tracy- Widom law. In addition, large deviation formulas for the 
maximum height are also computed. 



1 Introduction 



1.1 Background results 

It is a well-known result that U(iV) lattice QCD in two dimensions with Wilson's 
action [51] exhibits a third order phase transition in the large N limit [231 E3] . This 
is shown by forming the partition function for the plaquettes, which factorizes as a 
product of partition functions for each individual plaquette. The latter is identified 
with a zero-dimensional unitary matrix model having partition function given by 

G N (b) := (e bN ^ ut A , (1) 

where the matrices U G U(iV) are chosen with Haar measure and b is the scaled 
coupling. 

The matrix integral ([I]) depends only on the N eigenvalues of U, and in terms 
of these variables it can be written 

G N (b) = / de,... de N \{ e 2bN ™°> n l c< * - I 2 - ( 2 ) 

This can be interpreted as a partition function for a classical gas of charged par- 
ticles, confined to the unit circle, and repelling via logarithmic pair potential 
-(l/2)log|e ie - e**\ at the inverse temperature f3 — 2. The charges are also 
subject to the extensive one-body potential bN cos 9. In the form ^ the iV — > oo 
limit can be computed with the result 



iTooi^ logGjv(6) = { 26-|-|log26, 6>|, <5 (3) 

which is indeed discontinuous in the third derivative at b = 1/2. 

Some fifteen years after the work [231 E3] the same matrix integral Q appeared 
in a completely different setting. Consider a unit square, and place points uni- 
formly at random, with the number of points n chosen according to the Poisson 
distribution P(n) = ^j-e~ x2 with mean (n) = A 2 . Starting at (0,0) and finishing 
at (1, 1) form a piecewise linear path by joining dots with line segments of positive 
slope. There are evidently many such paths. From them, choose a longest path, 
i.e. a path with the maximum number of dots on it (see Fig. Q). Let h n denote 
the length of this longest path. Clearly hP is a random variable that fluctuates 
from one configuration of points to another and its probability distribution has 
been studied extensively in the context of the directed last passage percolation 
model due to Hammersley (see e.g. [HI §10.9]). For the cumulative distribution 
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of h n , one gets [2"Tl |4"5] (the first of the references gives a Toeplitz integral form 
equivalent to ([2]), while the second identifies the matrix integral explicitly) 

Pi(h D <N) = e- x2 (e XT < u+W A . (4) 

\ / ue\J(N) 

In the limit of large A, we set A = bN, it follows from ^ that in the large N limit 

Pr(/i D < N) — | e V 2 (-b 2 +2fe-3/4-(l/2)log26) £ > 1 ^ (5) 

hence providing for (|3| the interpretation as a large deviation formula for a prob- 
abilistic quantity in a statistical model [27] . 



(1, 1) 




(0, 0) 

Figure 1: Points (black dots) are distributed uniformly in a unit square with mean 
density A 2 . The length of any up-right path connecting the origin (0, 0) and the 
corner (1, 1) is defined by the number of points on the line. A path with the longest 
length is shown by the (blue) solid line. 



The desire to relate 2d lattice QCD to string theory focussed attention on the 
so-called double scaling limit of matrix integrals. Here, in addition to iV — > oo 
the coupling is tuned in the neighbourhood of the critical point b = 1/2 to give 
a well-defined scaling limit. It turns out that if one zooms in the neighbourhood 
of the critical point 6=1/2 and magnifies it by a factor iV 2 / 3 , i.e., one takes the 
limit (1/2 - b) 0, N -»■ oo, but keeping the product t = 2 4 / 3 (l/2 - b)N 2 ^ 3 fixed, 
then the second derivative of the free energy (the specific heat) tends to a function 
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of the single scaled variable t. In the case of the matrix integral (flj this double 
scaling limit was first analyzed by Periwal and Shevitz [12] whose result, using our 
notations, can be translated in the following form: setting 6 = 1/2 — 2 _4 / 3 iV~ 2 / 3 1 

^ lim [log G N (b) - b 2 N 2 ] = -g 2 (t), (6) 

where q(t) satisfies the special case a = of the Painleve II differential equation 

u" = 2u 3 + tu + a. (7) 

In other words 

q"(t) = 2q\t)+tq(t). (8) 

However no boundary condition was specified until Gross and Matytsin [22] refined 
the working of |42j to obtain a result which implies 

q(t) ~ Ai(t) (9) 

t— >oo 

(see also the earlier reference [39] for identification of Ai(t) in a similar context, 



and Section 3.2 below), where Ai(t) denotes the Airy function. 

In the context of the Hammersley model, recalling that in Q A = bN, we see 
that as a consequence of ^ 

lim Pr( -ft- ~ < t] = exp ( - [ (s - t)q 2 (s)ds) . (10) 



N^oc \(iV/2) 1 /3 

This result was first obtained in the probability literature [3] independent of the 
working of [12]. In fact, the authors of [3] were interested in studying a 'de- 
Poissonized' version of the Hammersley model, where the number of dots in the 
unit square is a fixed number N, and not a Poisson distributed random variable. 
This 'de-Poissonized' Hammersley model is, in turn, related to the so-called Ulam 
problem where one studies the statistics of the length of the longest increasing 
subsequence of a random permutation of iV integers (for a survey see [H El])- 
The length of the longest path in the 'de-Poissonized' Hammersley model has the 
same probability distribution as the length of the longest increasing subsequence in 
the random permutation. The Ulam problem has recently turned out to be a key 
model that connects various problems in combinatorics, physics and probability [U 



Remarkably, the right-hand side (rhs) of (10) admits a second interpretation 
within random matrix theory. Consider the Gaussian unitary ensemble (GUE) of 
the set of N x N complex Hermitian matrices X with measure proportional to 
e~ TrX . Let A max denote the largest eigenvalue. Its average, in the large iV limit, 
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is simply (A max ) ~ y2N. The typical fluctuations of A max around its average are 
very small of order ~ iV -1 / 6 . It turns out that the probability distribution of these 
typical fluctuations have a limiting distribution j4"9] 



hn^Prf 2 1 / 2 JV 1 / 6 (A max - V2N) < tj = exp ( - jf" (a - t)g 2 (s)ds^ := F 2 (t), 

(11) 

known as the /3 = 2 Tracy- Widom distribution. This distribution function has 
recently appeared in a number of problems ranging from physics to biology [4T1 134"] . 



So (10) and (11) link distribution functions in two seemingly unrelated problems. 
In addition, the same Tracy- Widom function also appears in the scaled specific 
heat of the U(iV) lattice QCD in two dimensions as demonstrated by Eq. (Jsl) . 



1.2 Statement of problems and summary of new results 

In the previous subsection, we have seen that the partition function of a two- 
dimensional field theory model [U(iV) lattice QCD with Wilson action], when mul- 
tiplied by a factor e~ A (see Eq. (4)), can be interpreted as the cumulative probabil- 
ity distribution of a certain random variable in the statistical physics/probability 
problem of Hammersley's directed last percolation model. It is then a natural 
question to ask if such connections can be established between other field theory 
models (on one side) and statistical physics models (on the other side). In this 
paper we establish another connection, namely between the continuum Yang-Mills 
gauge theory in two dimensions on a sphere (the field theory model) and the system 
of N non-intersecting Brownian motions on a line segment [0, L] (the statistical 
physics model). 

Non-intersecting random walkers, first introduced by de Gennes [14], followed 
by Fisher [18], have been studied extensively in statistical physics as they appear 
in a variety of physical contexts ranging from wetting and melting all the way to 
polymers and vortex lines in superconductors. Lattice versions of such walkers have 
also beautiful combinatorial properties [32]. Non-intersecting Brownian motions, 
defined in continuous space and time, have also recently appeared in a number of 
contexts, in particular its connection to the random matrix theory have been noted 
in a variety of situations [T7J [2E1 I2H HH 1131 El EI]- In this paper we introduce 
three new models of non-intersecting Brownian motions and establish their close 
connections to the Yang-Mills gauge theory in two dimensions on a sphere. 

Specifically, we consider a set of N non-intersecting Brownian motions on a 
finite segment [0, L] of the real line with different boundary conditions. Assuming 
that all the walkers start from the vicinity of the origin, we then define the reunion 
probability as the probability that the walkers reunite at the origin after a fixed 
interval of time which can be set to unity without any loss of generality. Next 
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we 'normalize' this reunion probability in a precise way to be defined shortly 
In one case, namely when both boundaries at and L are absorbing, one can 
relate this 'normalized' reunion probability to the probability distribution of the 
maximal height of N non-intersecting Brownian excursions. We show in this paper 
how to map this normalized reunion probability in the Brownian motion models 
to the exactly solvable partition function (up to a multiplicative factor) of two- 
dimensional Yang-Mills theory on a sphere. The boundary conditions at the edges 
and L select the gauge group of the associated Yang-Mills theory. We consider 
three different boundary conditions: periodic (model I), absorbing (model II) and 
reflecting (model III) which correspond respectively to the following gauge groups 
in the Yang-Mills theory: (I) periodic — > U(iV) (II) absorbing — > Sp(2iV) and (III) 
reflecting SO(2iV). 

Using the known results on the partition function from the Yang-Mills theory, 
we will show how these normalized reunion probabilities in the Brownian motion 
models can be related to the limiting Tracy- Widom distribution of the largest 
eigenvalue in some particular random matrix ensembles. The latter, in addition 
to the distribution J-^ify relating to complex Hermitian matrices, involves its com- 
panion T\(t) for real symmetric matrices. Explicitly with the GOE specified as the 
set of N x N real symmetric matrices X with measure proportional to e~ TrX / 2 , 
and A max denoting the largest eigenvalue, one has [50] 



Jim^Pr(V2iVi(A max - V2N) < ?j = exp ^ - ((s - t) q 2 {s) - q(s)) dsj 

■= ?i{t). (12) 

We consider the following three models of N non-intersecting Brownian walkers 
on a one- dimensional line segment [0, L] with different boundary conditions. Let 
the walkers be labelled by their positions at time r, i.e., by Xi(r) < ... < xn(t). 

Model I: In the first model we consider periodic boundary conditions on the 
line segment [0, L\. Alternatively, one can think of the domain as a circle of 
circumference L (of radius L/2tt). All walkers start initially in the vicinity of a 
point on the circle which we call the origin. We can label the positions of the 
walkers by their angles {6 1 ,6 2 , ■ ■ ■ ,6^} (see Appendix A for details). Let the 
initial angles be denoted by {ei,e2, . . . ,ev} where e^'s are small. Eventually we 
will take the limit — > 0. We denote by -R£(l) the reunion probability after time 
t = 1 (note that the walkers, in a bunch, may wind the circle multiple times), 
i.e, the probability that the walkers return to their initial positions after time 
t = 1 (staying non-intersecting over the time interval r G [0, 1]). The superscript 
/ corresponds to model I. Evidently i?£(l) depends on N and also on the starting 
angles {ei, €2, ■ ■ ■ , cat}- To avoid this additional dependence on the ej's, let us 
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introduce the normalized reunion probability defined as the ratio 



G N {L) = , (13) 



oo v 



where we assume that we have taken the 6j — > limit. The e dependence actually 
cancels out between the numerator and the denominator (see Appendix A) and 
hence GV(L) depends only on iV and L. In Appendix A, we calculate GV(L) 
explicitly and show that 

.00 00 

Gn{L) = ^ E • • E AK-.n^^^, (14) 

ni=— 00 njv=— 00 

where 

A(n h n N ) = ] [ (rii - rij) , 



l<i<j<N 



and the prefactor 



A N = = AT — 1 (15) 



1 

(2n)^ 2 T(j + 2) 

ensures that Gn{L — > 00) = 1. In the next section we will see that this normalized 
reunion probability Gn(L) is, up to a prefactor, exactly identical to the partition 
function of the 2-d Yang-Mills theory on a sphere with gauge group U(iV). 

We remark that for the non-intersecting Brownian motions on a circle, a similar 
mapping was first noticed by Minahan and Polychronakos [38] , with a slightly dif- 
ferent normalization than ours. However, the behavior of the normalized reunion 
probability G^{L) as a function of the system size L was not analysed in [38] and 
consequently they did not uncover the existence of the Tracy- Widom distribution 
^(t) near the critical point L C (N) = 2\/N for large N in the reunion probability, 
which is indeed one of our main findings in this paper. 

Model II: In the second model the domain is the line segment [0, L] with ab- 
sorbing boundary conditions at both boundaries and L. Once again, the N non- 
intersecting Brownian motions start initially at the positions, say, {ei, e 2 , . . . , e^} 
in the vicinity of the origin where eventually we will take the limit q — > for 
all i. As in Model I, we define the reunion probability -R£ 7 (l) as the probability 
that the walkers return to their initial positions after a fixed time r = 1 (staying 
non-intersecting over the time interval r G [0, 1]). Analogous to Model I, we define 
the normalized reunion probability 

F,(L) = g|) , (16) 
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which becomes independent of the starting positions q's in the limit when e« — > 
for all i. Hence, F^(L) depends only on N and L. 

This ratio Fn{L) in Model II also has a different probabilistic interpretation. 
Consider the same model but now on the semi-infinite line [0, oo] with still absorb- 
ing boundary condition at 0. The walkers, as usual, start in the vicinity of the 
origin and are conditioned to return to the origin exactly at r = 1 (see Fig. ([2])). 
If one plots the space-time trajectories of the walkers, a typical configuration looks 
like half of a watermelon (see Fig. ([2])), or a watermelon in presence of an absorb- 
ing wall. Such configurations of Brownian motions are known as non-intersecting 
Brownian excursions and their statistical properties have been studied quite exten- 
sively in the recent past. A particular observable that has generated some recent 
interests is the so-called 'height' of the watermelon [6J HU EDI USB HZl EI] defined as 
follows (see also Ref. [8] for a related quantity in the context of Dyson's Brownian 
motion). Let denote the maximal displacement of the rightmost walker xn 
in this time interval r G [0,1], i.e., the maximal height of the topmost path in 
half-watermelon configuration (see Fig. ([2])), i.e., H N = max r {xAr(r), < r < 1}. 
This global maximal height is a random variable which fluctuates from one 
configuration of half-watermelon to another. What is the probability distribution 
of H^l For N = 1 the distribution of Hn is easy to compute and already for 
iV = 2 it is somewhat complicated |30| . Recently, however, an exact formula for 
the distribution of H^, valid for all N, was derived in [UJ using Fermionic path 
integral method. 

To relate the distribution of in the semi-infinite system defined above to the 
ratio of reunion probabilities in the finite segment [0, L] defined in (16), it is useful 
to consider the cumulative probability Pt(Hn < L) in the semi-infinite geometry, 
where L now is just a variable. To compute this cumulative probability, we need 
to calculate the fraction of half-watermelon configurations (out of all possible half- 
watermelon configurations) that never cross the level L, i.e., whose heights stay 
below L over the time interval r G [0,1] (see Fig. ([2])). This fraction can be 
computed by putting an absorbing boundary at L (thus killing all configurations 
that touch/cross the level L). It is then clear that Pt(Hn < L) is nothing but the 
normalized reunion probability F^{L) defined in (16). As mentioned above, this 
cumulative probability distribution of the maximal height was computed exactly 
in Ref. [U] (see also [IE], [31] for related computations) 

F N (L) := Pr(H N < L) 

„ oo oo N 

AV,.../<)(n<K &Ef -"' • (in 



l / 2N 2 +N 

ni=l rijv=l 



where B N = ix 2N2 + N / (2 n2 - n I 2 \\^ T{2 + j)T(3/2 + j)). A brief derivation of 
this result is given in Appendix B. Note that a similar probabilistic interpretation 
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(i.e. as a cumulative distribution of a random variable) does not exist for the 
normalized reunion probability Gn(L) in Model I. This is evident from the fact 
that, unlike -Fjv(L) in Model II, the quantity Gn(L) in Model I is not bounded 
from above by 1 (this can even be seen by a direct computation in the N = 1 case, 
see (102) in Appendix A). 




Figure 2: Trajectories of N non-intersecting Brownian motions X\(r) < ^(t) < 
... < Xatt), all start at the origin and return to the origin at r = 1, staying 
positive in between. F^(M) denotes the probability that the maximal height 
Hn = max T {xAr(r), < r < 1} stays below the level L over the time interval 
< r < 1. 



In the next section, we will see that the normalized reunion probability (or 
equivalently the maximal height distribution) Fn(L) in Model II is, again up to 
a prefactor, precisely equal to the partition function of the 2-d Yang-Mills theory 
on a sphere, but with a gauge group Sp(2iV) (as opposed to the group U(iV) in 
Model I). 

Model III: We consider a third model of non-intersecting Brownian motions where 
the walkers move again on a finite line segment [0, L], but this time with reflecting 
boundary conditions at both boundaries and L. Again the walkers start in the 
vicinity of the origin at time r = and we consider the reunion probability i?^ /7 (l) 
that they reunite at time r = 1 at the origin. Following Models I and II, we define 
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the normalized reunion probability 

E N (L) = jgg , (18) 

that is independent of the starting positions {ei, €2, ■ ■ ■ , e^v} in the limit when all 
the ej's tend to zero and hence depends only on N and L. Following similar steps 
as in Models I and II, but with reflecting boundary conditions at both and L, 
we find the exact expression (see Appendix B) 



En(L) = j^ £ - £ AV,...,n? r ) e -^)^»5 ) (19) 

m=— oo njv=— 00 



£2V 2 -iV 

and the prefactor 



^2N 2 -N 23N/2-N 2 

n?^r(2 + j)r(i/2 + j) 



^ = ^^ivri^-^TTT^-V (20) 



ensures that E N (L — > 00) = 1. In the next section we will see that E N (L), up 
to a prefactor, is exactly identical to the partition function of the 2-d Yang-Mills 
theory on a sphere with gauge group SO(2iV). 

Thus, we find that by changing the boundary conditions at the edges of the 
line segment [0, L] in the non-intersecting Brownian motion models we can relate 
the normalized reunion probability to the partition function of the 2-d Yang-Mills 
theory on a sphere with an appropriate gauge group. Model I, II and III corre- 
spond respectively to gauge groups U(iV), Sp(2iV) and SO(27V). 

Summary of new results: Let us then summarize the main new results in this 
paper: 

• We have shown that the normalized reunion probability of a set of N non- 
intersecting Brownian motions moving on a line segment [0, L\ with a prescribed 
boundary condition is identical, up to a prefactor, to the partition function of the 
2-d Yang-Mills theory on a sphere with an appropriate gauge group which depends 
on the boundary conditions in the Brownian motion model. We have shown that 
three different boundary conditions lead respectively to the gauge groups: (I) pe- 
riodic ->■ U(N) (II) absorbing ->■ Sp(2iV) and (III) reflecting ->• SO(2iV). 

• The partition function of the 2-d Yang-Mills theory with a given gauge group 
is exactly solvable and many results are known. In particular, it is known that in 
the so-called 'double scaling' limit, the singular part of the Yang-Mills free energy 
satisfies a Painleve equation. To use these results for the Brownian motion model 
via the correspondence established above, we need to however treat the prefactor 
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of the correspondence properly in the double scaling limit. Taking into account 
new terms arising via these prefactors, we show that the normalized reunion prob- 
abilities, appropriately centered and scaled as a function L for large but fixed N, 
also share a 'double scaling' limit where they are precisely described by the Tracy- 
Widom distribution. In the periodic case, one gets ^(t) while for the absorbing 
and reflecting cases one gets J r i(t). Thus, this relates for the first time (to our 
knowledge) the Painleve equation that appears in the Yang-Mills free energy to 
the Painleve equation in the Tracy- Widom distribution. 

• As a byproduct of this mapping, we also show that the 3-rd order phase tran- 
sition between the strong and the weak coupling phases (separated by the double 
scaling regime) in the 2-d Yang-Mills theory translates into a 3-rd order phase 
transition in the behavior of the normalized reunion probability in the Brownian 
motion model, as one varies the system size L across a critical value L C (N) ~ \/N 
for large but fixed N. The strong and weak coupling phases correspond respec- 
tively to the left and right large deviation (away from the critical value L C (N)) 
behaviors of the normalized reunion probability as a function of L. Again using 
results from the Yang-Mills theory (taking into account the prefactors correctly), 
we thus obtain precise large deviation behaviors of these reunion probabilities. In 
particular, in the right large deviation behavior, we find a new type of crossover 
phenomenon. 

• For the special case of absorbing boundary condition (Model II), this gives 
us a direct proof that the distribution of the maximal height Hn of a set of iV 
non-intersecting Brownian excursions, when properly centered and scaled, has the 



Tracy- Widom distribution JF\(t) described in (12). This was shown before rather 
indirectly in Ref. [2H] via a mapping to a polynuclear growth model. Here we 
obtain a direct proof of this result. 

Let us remark that while some aspects of the analogies between non-intersecting 
Brownian paths and Yang-Mills theory on the sphere have been noticed in earlier 
publications [26, 25J, this precise correspondence between the normalized reunion 
probability in the Brownian motion models (with different boundary conditions) 
and the partition function of the 2-d Yang-Mills theory on a sphere (with different 
gauge groups) seems to be new, to the best of our knowledge (except for the peri- 
odic case when a similar correspondence was noted in |38j). More importantly, the 
probabilistic connection between the Yang-Mills partition function in the double 
scaling limit and the Tracy-Widom distribution of the largest eigenvalue of a ran- 
dom matrix (established here using the connection via non-intersecting Brownian 
motions), to our knowledge, was not noticed earlier. 
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The rest of the paper is organised as follows. In Section 2, we briefly recapitu- 
late the exact solution of the continuum Yang- Mills theory in two dimensions and 
then establish the correspondence between the partition functions of the gauge the- 
ory with normalized reunion probabilities in the Brownian motion models defined 
above. Next we study the consequences of this correspondence for Model I and 
Model II respectively in Sections 3 and 4. In particular, we will see how the Tracy- 
Widom distributions J-zit) and F\(t) emerge in the double scaling limit in Models 
I and II respectively. For Model II, this correspondence also provides us with 
the precise asymptotic distribution of the maximal height for iV non-intersecting 
Brownian excursions. We do not discuss Model III in details here as the analysis 
is very similar to that of Model II. The detailed derivation of the expression of 
Gn(L) in (14) is provided in the Appendix A. The derivations of F N (L) in (17) 
and of En(L) in (19) follow via similar calculations which are briefly outlined in 
Appendix B. 



2 Correspondence between 2-d Yang-Mills the- 
ory and non-intersecting Brownian motions on 
a line segment [0, L] 

We start by briefly recapitulating how one computes the 2-d Yang-Mills partition 
function [23 SSI 121] (see also the review by Cordes, Moore and Ramgoolam 0). 
Consider a general orientable two-dimensional manifold M. with corresponding 
volume form yjg. At each point x of this manifold sits a gauge field A^(x) (with 
the space index fi = 1, 2) which is an (N x N) matrix. For simplicity, let us first 
consider pure U(iV) Yang-Mills gauge theory whose partition function in contin- 
uum is defined by the functional integral 

Z M = J [vA^e-W^lMVsn^F^ ) (21) 

where A is a coupling constant and the field strengths are defined by 

= d^Ay - d v A^ + iL4 M , A v \ . (22) 



Under a local gauge transformation A^ — > S~ 1 (x)A fM S(x) — iS~ 1 (x)d^ 1 S(x) (where 
S(x) is an (N x N) unitary matrix), the field strengths transform as F^ — > 
S~ 1 (x)F flu S(x) thus keeping the action in (21) gauge invariant. 

The partition function in (21) in two dimensions can be computed exactly via 

S3. 



the original idea due to Migdal [37] . One can actually use a particular lattice 
regularization of the continuum theory which is both exact and additive in the 
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following sense [24] . One can divide the manifold into polygons (for example one 
can choose triangles as basic units or plaquettes) and define a unitary matrix Ul 
sitting at the center of each link L of this triangulated manifold (see Fig. [3]). Then 
the lattice regularized partition function can be written [23] 

(lattice) = J H dU L J] Zp[U P ], (23) 

L plaquettes 

where Up = riiepiaquctte Ul is called the loop product introduced by Wilson [53] 
and Zp{Up) is some appropriate lattice action associated with the plaquette P. 
The only constraint on the choice of Zp is that it should reduce to the continuum 
action when the plaquette size goes to zero. 

A standard choice for Zp(Up) is the Wilson action 



Zp(U P ) = exp \bNTr(U P + U P )] , (24) 
which does reduce to the continuum action in the limit when the plaquette area 



goes to zero. With this choice, the lattice partition function in (23) is exactly 
solvable [231 [53] as it reduces to computing a single matrix integral in Eq. (JT1) 
already discussed in the introduction. However, the Wilson action is not invariant 
under renormalization in the following sense. Following Migdal [37], one can take 
two adjacent plaquettes Pi and P 2 with their respective actions Zp 1 and Zp 2 and 
fuses them to form a bigger plaquette with area equal to that of Pi + P 2 , after 
integrating out the unitary matrix sitting on the common link between the two 
plaquettes (see Fig. [3]). This gives the Migdal recursion relation 

J dU 3 Z Pl (U 1 U 2 U 3 )Z Fil (U 4 Usl4) = Z' Pi+P2 {U x U 2 U A U 5 ) . (25) 



In general, the renormalized plaquette action Z' p does not have the same func- 
tional form as the bare action Zp (which is indeed the case when one chooses 
Wilson action as the bare action). So, the natural thing to look for is the fixed 
point solution of this recursion relation that keeps the form of Zp invariant un- 
der renormalization. Indeed, Migdal [37] and later Rusakov [3S] showed that the 
appropriate fixed point action is given by the so-called heat-kernel action 



Z P = ^2d R x R (U P )exp 



R 



XAj 



2N 



C 2 (R) 



(26) 



where the sum runs over all irreducible representations R of U(iV), dp is the 
dimension of R, Xr(Up) is the character of Up in the representation R, C 2 (R) is 
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Figure 3: a typical triangulation of the two-dimensional manifold with unitary 
matrices Uls on the edges. One can fuse two triangles P\ and Pi by integrating 
the matrix U3 along their common edge and obtain a parallelogram with four edges 
with matrices Ui, U 2 , U4 and U 5 on these edges. 



the quadratic Casimir operator of R and Ap is the area of the plaquette. The 
coupling constant A is fixed and will henceforth be chosen to be unity: A = 1. One 
can also verify that this heat-kernel action reduces to the continuum action in the 
limit when the area of each basic plaquette goes to zero [37J SB] • Note that this fixed 
point choice of Zp makes the lattice representation of the continuum theory exact 
in the sense that the final result is independent of the triangulation as one can add 
as many triangles (and in whichever way) to cover the manifold thus approaching 
the continuum limit [24J. The name heat-kernel stems from the fact that one can 
express the heat-kernel action Zp(U) = (U\ exp[— (l/2iV)AA]|l) where A is the 
Laplacian on the group [36J . This fact already gives a hint that there might be an 
underlying diffusion process inbuilt in this effective action, though the precise fact 
that they correspond to non-intersecting Brownian motions is still not evident at 



this point. Using the heat-kernel action Zp in (26) one can then evaluate the full 
partition function on the manifold [46J 



Z M = S ^d 2 R 2v exp 



ft 



(27) 



where p is the genus of the manifold and A is the total surface area. This then 
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represents an exact solution of the partition function of the continuum Yang-Mills 
theory on a twodimensional manifold. Note that even though we have specifically 



used the group U(iV) for the discussion above, the result (27) is valid as well for 
other groups such as SU(iV), Sp(iV) and SO(iV). 
For a sphere, using p = 0, one gets 



Zm = S ^d\ ex P 

R 



A 
~2N 



C 2 (R) 



(2f 



The irreducible representations can be labelled by the lengths of the Young 
diagrams and one can explicitly express the partition function as a multiple sum. 
For the group U(iV), the partition function reads 



Z M :=Z{A;U(N)) = c N e 



2 °° 
vT -1 A 9/ \ 

24 2^ A z [ni,...,n N )e 

rii , . . . , n — — oo 



(A/2iV)Ef =1 ™ 2 



(29) 

where A(n%, ri2, ■ ■ ■ , tin) is the van der Monde determinant defined in ( 1.2[ ) and 
is a constant independent of A. For the groups Sp(2iV) and SO(2iV) one can 
similarly express the partition function as a multiple sum [11]. For example, for 
Sp(2iV) one gets 



oo N 

Z(A; Sp(2iV)) = c N e A ^ N+1 ^ £ A 2 « ...,<)( J] n^e~^ 



ni,...,njv=— oo 



(30) 

where cn is independent of A. Similarly, for the SO(2iV) group, one obtains [TT] 

oo 

Z(A;SO(2N)) = b N e A ^ N -^ ^ A 2 (n?, . . . , n 2 N )e~^ n " , (31) 



ni,...,rajv— — oo 



where 6 at is independent of A. 
Comparing the formulas in 



14), (17) and (19) with the expression for partition 



functions respectively in (29), (30), and (31), we see that the normalized reunion 



probabilities in the three models of non-intersecting Brownian motions, up to pref- 
actors that can be computed explicitly, are isomorphic to the partition functions 
of the Yang-Mills theory on a sphere with respective gauge groups U(iV), Sp(2iV) 
and SO(2iV), provided we make the identification A/2N 2n 2 /L 2 in Model I, 
A/ AN — > 7r 2 /2L 2 in model II and A/ AN — > ix 2 / L 2 in model III. More precisely, up 
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to known prefactors, 



G N (L) oc Z[A=— ;U(JV)J , (32) 



F N (L) oc Z(A= 2 ^;S V (2N)) , (33) 



/ 4n 2 N \ 
E N (L) oc ZU = -^;SO(2iV)J . (34) 

This correspondence between the normalized reunion probability for non-intersecting 
Brownian motions with different boundary conditions and the Yang-Mills parti- 
tion functions on a sphere with corresponding gauge groups is one of the main 
observations of this paper. 

In the next two sections we study the consequences of this correspondence in 
detail for Model I and II. We skip detailed studies of Model III since it can be 
handled exactly in the same way as Model II. The main point is to derive the 
precise asymptotics of the normalized reunion probabilities in the two models (in 
particular for Model II this will give us the asymptotic behavior of the distribution 
of the maximal height for non-intersecting Brownian excursions) by using the 
known behavior of the asymptotic properties of the partition functions of the 
corresponding gauge theory, albeit taking into account correctly the L and N 
dependence of the respective prefactors. 



3 Brownian walkers on a circle: Model I 



Comparing (29) and (14) we have the following exact identity between the normal- 



ized reunion probability and the U(iV) Yang-Mills partition function on a sphere 
(upon substituting A = 4ir 2 N/L 2 in p9~])) 



G N (L) 



A 



N e -K 2 N(N 2 -l)/6L 2 j^-N 2 
CN 



Z 



4vr 2 iV 



TJ(A0j 



(35) 



where the constants A^, given in (15), is independent of L. Similarly the constant 



is independent of L and can be fixed in any way. Later, we choose cn such that 
for large N, log cn — —N 2 logN. This choice of cn ensures that the free energy 
\ogZ M ~ 0{N 2 ) for large N (see below). 

In Ref. [U] it was shown that for large N, 2(A; U(iV)) exhibits a 3-rd order 
phase transition at the critical value A = n 2 separating a weak coupling regime 
for A < tx 2 and a strong coupling regime for A > ix 2 (see Fig. [4|. This is 
the so called Douglas-Kazakov phase transition |15j . Using the correspondence 
L 2 := 4n 2 N/A, we then find that Gn{L), considered as a function of L with N 
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G N {L) 



TW 



N 




right tail ^2 _ iw 2 N 



2VN 



A 
> 



L 



critical region 



weak 
coupling 



strong 
coupling 



'A 



7T 



Figure 4: Left: Schematic sketch of G^iL) as defined in Eq. (13) for iV vicious 
walkers on a circle, as a function of L, for fixed but large N. Right: Sketch 
of the phase diagram in the plane (A, N) of two-dimensional Yang-Mills theory 
on a sphere with the gauge group U(iV) as obtained in Ref. [15J. The weak 
(strong) coupling phase in the right panel corresponds to the right (left) tail of 
Gn(L) in the left panel. The critical region around A = n 2 in the right panel 
corresponds to the Tracy- Widom (TW) regime in the left panel around the critical 
point L C (N) = 2y/N. 



large but fixed, must also exhibit a 3-rd order phase transition at the critical value 
L C (N) = 2y/N. Furthermore, the weak coupling regime (A < n 2 ) corresponds to 
L > 2yN and thus describes the right tail of Gn(L), while the strong coupling 
regime corresponds to L < 2\/N and describes instead the left tail of Gn(L) (see 
Fig. [4|. The critical regime around A = it 2 is the so-called "double scaling limit" 
in the matrix model and has width of order iW 2 / 3 . It corresponds, as we will see, 
to the region (of width iW 1 / 6 around L = 2\^N) where Gn(L), correctly shifted 
and scaled, is described by the Tracy-Widom distribution ^(t). 



3.1 Large deviation tails 

Let us first consider the behavior of Gn(L), as a function of L, in the two tails 
far away from the critical value L C (N) = 2\/N for large but fixed N. To do 



this, we can exploit the exact identity in (35) and use the known properties of 
the U(iV) partition function respectively in the strong (A > it 2 ) and the weak 
(A < tc 2 ) coupling phases. The left (L < 2y/N) and the right (L > 2\fN) tails 
of Gn(L), away from the critical point L C (N) = 2v / W, correspond respectively to 
the behavior of the partition function in the strong (A > it 2 ) and weak (A < tt 2 ) 
coupling phases. 
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Let us begin by summarizing the known properties of the U(iV) partition func- 



tion in the strong and weak coupling phases. The summand in (29) can naturally 
be regarded as a function of rii/N (i = 1, . . . , N). In the large N limit the variables 
{ni/N}i = i r .. t N approximate the coordinates of a continuous N particle system. As- 
sociated with the particles is a density p(x), but because the lattice spacing is 1/N 
and there are N particles the density at any point cannot exceed 1, and thus 
p(x) < 1 for all x. By using this viewpoint to perform a constrained saddle point 



analysis of (29) (see also Section 4.1 below), it was shown by Douglas and Kazakov 
[15] that 

1 . _ f F-{A), A < 7T 2 , 



where 



fc^'oS^UW)-! F+(A), A>* 



Q A 1 

FJA) = -- log A. 

\ ) 4 24 2 s 



(36) 



(37) 



The function F + (A) about A = n 2 agrees with F_(A) up to 0((A — re 2 ) 2 ) in its 
power series expansion. Thus [T51 eq. (40)] 



F + {A)-F.{A) 



3tt' 



71 



2\3 



+ C((A-7T 



2^4^ 



(3f 



Its explicit form — or more precisely that of its derivative — is given in terms 
of elliptic integrals. Specifically, with K = K(k), E = E{k) denoting standard 
elliptic integrals, where k is specified by the requirement 



A 



(2E - k' 2 K)K , k' 2 = 1 - k 2 



(39) 



and 



one has [151 ec L- (35)] 



F'AA) 



a 



AK 

a 2 k' 2 _ 1 a A k' A A 
24 



12 



96 



(40) 



(41) 



In (36), the case A < tx corresponds to the weak coupling phase and the 



constraint on the density of the variables rii/N being less than 1 can be ignored, 
implying that the Riemann sum is well approximated by the corresponding multi- 
dimensional integral. In the latter the A dependence can be determined by scaling, 
and the resulting integral evaluated explicitly (see e.g. [TjJl eq. (4.140)]) to deduce 
the stated result. However for A > n 2 doing this would imply that the density is 
greater than 1, so the discrete sum is no longer well approximated by the continuous 
integral. The density saturates at 1 for some range of values of the variables rii/N 
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about the origin, but goes continuously to zero to be supported on a finite interval 
(see e.g. Fig. 2 in 0). 

Let us now use these known results on the partition function in our exact 
identity (35) to derive the corresponding large deviation properties of Gn{L). 
Setting, for convenience, L = 2y/Nr in (35) so that r = 1 corresponds to the 
critical point, we get 



A p vr 2 (Af 2 -l)/24r 2 . 2 



c N (2VNr) 



(42) 



We first note that A N in Eq. (Il5j contains the product r (i + 2 ) = G ( N + 2 ) 

in the denominator, where G(x) denotes the Barnes G-function whose asymptotic 
expansion is (see e.g., [HI eq.(4.184)]) 



ln(G(z + l)) 



z 2 In z 



3 4 + - 4 ln W + C'(-l) + 0(1/ z) . (43) 
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Using this result and the choice cn — e 
term cancels when one considers the ratio 

y 4 Are vr 2 (V 2 -l)/24r 2 



-N 2 log TV 



, we see that the leading N 2 log iV 



c N (2y/Nr) N2 



V 2 (3/4-log(r/ 7 r)+7r 2 /24r 2 )+C'(V) 



(44) 



Consequently, the result in Eq. (|36|) gives the large deviation formula 

G 7V (2v / iVr) r 



1, 



-N 2 {F-(7r 2 /r 2 )-F+(TT 2 /r 2 )) 



r > 1 
r < 1. 



(45) 



This is then a new result on the far tails of the normalized reunion probability 
Gn(L), as a function of L, for fixed but large N. 



Note that the precise meaning of ~ in (|45|) is the following 
\G N (2VNr) 



- lim ^jln 



0, r > 1 

F_(7r 2 /r 2 ) -F+(7r 2 /r 2 ), r < 1. 



(46) 



Thus this calculation only gives the leading 0(N 2 ) term. While in the left tail 
(r < 1) the leading term of (9(iV 2 ) is a finite nontrivial quantity, this leading 
0(N 2 ) term actually vanishes in the right tail (r > 1). So, to obtain the finer 
behavior in the right tail one needs to keep track of the next subleading correction 
of 0(N) term in \n[G N (2VNr)} for r > 1. 

In the gauge theory context, this means that we need to obtain in the weak 
coupling phase, not only the leading term of 0(N 2 ) but also the subleading cor- 
rections. Fortunately, these subleading corrections in the free energy lni^ in the 
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weak coupling phase were also calculated by Gross and Matytsin [2"2"j eq. (2.41) 



A (-1) N A ( A\~V* 2^N,As 



24 v^TriV 2vr 2 V vr 2 , 
+ O e -—^> , (47) 



up to terms independent of A (these will depend on the precise form of cn in (J29J 
where 

7 (x) = - £ log ] + V ^ = |(1 - a;) 3 / 2 + - a;) 5 / 2 ). (45 

z 1 — vl — ^ ^ 



We can then use this result in our identity (42) to predict the following right large 
deviation tail 

l-G N (2VNr) ~ (_i)^ e -2^7(i/r a ) j r>1? (49) 

which shows an interesting oscillatory behavior. Again this is a new result for 
the normalized reunion probability of the non-intersecting Brownian motions on 
a circle. Let us recall that Gn(L) does not have the meaning of a cumulative 
distribution and hence the oscillating sign of [1 — G^{2^/Nr)\ (with N) is not 
really problematic. 

3.2 Double scaling limit 

Having obtained the precise large N asymptotic behavior of G^{L) as a function 
of L in the left [L < 2\/N) and the right (L > 2yN) tails, we now turn our 
attention to the behavior of Gn{L) in the vicinity of the critical point, i.e., when 
L is close to L C (N) = 2y/N. In the gauge theory, this corresponds to the double 
scaling regime near the critical point A = tt 2 . Below, we first discuss the known 
results on the partition function in the double scaling regime. We then use these 



results in our exact identity (35) to show that correspondingly Gn(L), in a narrow 
region \L — 2y/N\ ~ jV -1 / 6 around the critical point L C (N) = 2\/N } has the scaling 
behavior 

G N (L) -> JS (2 2/3 N 1/6 \L - 2y/N\) , (50) 

where the scaling function ^(t) is precisely the Tracy- Widom distribution function 
for GUE random matrices defined by 

? 2 {t) = exp [ - [ (s- t)q 2 (s)ds] , (51) 
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with q(s) satisfying the Painleve equation (j8j). This is the main result of this 
subsection and details are provided below. 

Gross and Matytsin [22] used the method of orthogonal polynomials (see e.g. [191 



Ch. 5]) to analyse the N — > oo asymptotic behaviour of (29) when the coupling A 
is tuned in the neighbourhood of the transition point A = tt 2 . Explicitly, in what 
is referred to as the double scaling limit, this is achieved by taking iV — > oo while 
keeping (tt 2 — A)N 2 ^ fixed. 

As most clearly set out in [TU eq. (23)], for certain orthogonal polynomial 
normalizations , the method of orthogonal polynomials gives 

^logZ(A;U(iV)) = ^R^iR^l + R^) (52) 

with + (— ) chosen according to N being even (odd). Moreover, it is shown in [22] 
that for j near N 



AT 



2 



R ? ] = ~2 ^ (-l)^ 5/3 /i(^) + 0(iV 4/3 ) (53) 

7T 

where 

X = N 2 '\l-^). (54) 



In (53) fi satisfies the differential equation 



f" ~ 4xfi ~ y /i 3 = , (55) 



subject to the boundary condition 



/ 2 e~i x 1 

fx(x) ~ -y— — , (56) 



(here we have corrected a factor of 7r 2 in Eq. (5.13) of [22] in writing (55), while 



in writing (56) we have changed the sign in Eq. (5.14); to see that the latter is 
needed compare Eq. (5.14) with Eq. (5.4)). 

As noted in [22], (55) can be identified with the Painleve II equation ^ in the 



x = 2~ 2 l\ f l ( x ) = -— u (t) (57) 



case a = 0. Explicitly, the latter is obtained upon the substitutions 

2 5/3 
TT 

and furthermore, upon recalling the asymptotic expansion 

e -2t 3 / 2 /3 

Ai(t) tZo 27^ ' 
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a boundary condition consistent with ^ is obtained. 

Substituting (53) in (52) implies that with x as specified by (54) fixed 
below (5.16)] 

d 2 N 2 / 2r 7r 4 

logZ(A;U(iV)) ~ — fl - _ - — ^A 2 (x) 



dA 2 ta v ' v " 2vr 4 V N 2 / 3 2N 2 / 3 



or equivalently, upon recalling (54) and making use of the variables (57), 

^\oge- A2N2 ^ 4 Z(A;XJ(N)) = -q 2 (t) , (58) 

where q(t) is as i n (|6~| ). 

The relation (|42[) now tells us that for the ratio of return probabilities we have 
that for iV — > oo 

^ log (5^(2^(1 + t/2(2iV) 2/3 )) = -q 2 {t). (59) 



But the distribution (11) for the scaled largest eigenvalue in the GUE satisfies this 
same relation, and so we have 

lim G N (2yfN{l + t/2(2iV) 2/3 )) = F 2 (t). (60) 

N— too 



Here use has also been made of the fact, which follows from (45), that the left-hand 
side (lhs) tends to 1 as t — > oo. 

We know from previous studies of large deviation formulas associated with the 
largest eigenvalue of a random matrix [T21 H31 E21 ESI EH Ej, where the transition 
region is specified by the Tracy- Widom scaling function, that the expansion of 
the large deviation functions around the transition point coincides with the tail 
behaviours of the transition region Tracy- Widom scaling function. Interestingly, 
this property holds in the present setting for one side (left) of the tails only. Thus 
making use of the expansions ( [38] ) and ( [48] ) we obtain 

a 

t — y — oo 



G n (2VN(1 + t/2(2iV) 2/3 ) ~ eT5 (61) 



G' n (2VN(1 + t/2(2iV) 2/3 ) ~ (62) 



t— >oo 



with the first of these only the precise tail form J-^if) [IH]. The t — > oo tail of 

_4t£/2 

has the leading form e 3 , and so differs by a factor of 2 in the exponent. This is 
the indication of an interesting crossover which we describe here qualitatively |55j . 
Close to L = 2VN, with L - 2^/N ~ 0(N~ 1 ^ 6 ), a calculation of G N (2y/N(l + 
t/2(2iV) 2 / 3 )) beyond leading order shows that 

hgG N (2y/N(l + t/2{2N) 2 ^)) = log (T 2 (t)) + N^-lf^t) + o^ 1 / 3 ) ,(63) 
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where the function g(t) can be expressed explicitly in terms of q(t) and behaves 
asymptotically as 



git) ~ e 3* 



when £ — 7- oo 



(64) 



What happens when one increases L from the critical region L — 2\/N ~ (^(iV -1 / 6 ) 
towards the large deviation regime in the right tail, L > 2y/N (see Fig. ^ ? As L 
increases away from L C (N) = 2\/N, the amplitude of the second term in the rhs 
of Eq. (63), which is oscillating with N, increases relatively to the amplitude of 
the first term. And at some value L = L cmss (N), it becomes larger than the first 
one: in the large deviation regime it becomes the leading term (still oscillating 



with N), as given in Eq. (49). On the other hand, the first term in the rhs of Eq. 



(63), log (J-2(t)), is subdominant for L > L cvoss (N) and becomes actually the term 

_ 4w 2 N ( A \ 

of order 0[e A 1 \^ z >) in Eq. (47). Hence, there is a crossover between the two 
terms in the rhs of Eq. (63) as L crosses the value L cvoss (N). Balancing these two 



terms and making use of the leading behavior of the right tail of JFzit) together 



with the asymptotic behavior of g(t) given in Eq. (64), one obtains an estimate of 



L CTOSS (N) as 



,(JV) - 2\fN ~ N- 1/6 (\ogN) 2/3 . 



(65) 



Note that such a crossover is absent in the distribution of the largest eigenvalue 
of GUE random matrices and it is thus a specific feature of this vicious walkers 
problem. 



4 Brownian walkers with absorbing walls: Model 
II 

In Model II, we have non-intersecting Brownian motions on [0, L] with absorbing 



boundary conditions. Comparing (17) and (30) we have the exact identity 



F N (L) = ^ e ^+V2)(iV + i))/6^ l .2N--n z ( to?N Sp{2N) \ (66) 
c N \ L 2 J 

where the constants Bn and &n are independent of L. As in the U(iV) case, 
thanks to this identity and known large N properties of the partition function 
Z (A; Sp(2iV)) provides us new large N results for normalized reunion probability 
F N (L) as a function L. In addition, in this case, F N (L) is also identical to the 
distribution of the maximal height for non-intersecting Brownian excursions. 
We thus get, as a bonus, new asymptotic results for the height distribution: both 
in the critical regime where its scaling behaviour is described by the Tracy- Widom 



23 



distribution T\(t) of GOE matrices, as well as in its large deviation tails. As in 
the previous section, below we first discuss the large deviation tails and then the 
behavior in the critical regime. 



4.1 Large deviation tails 

In this case, using the correspondence A = 2n 2 N/L 2 , the critical point A = it 2 in 
the gauge theory corresponds to a critical value L = y/2N . For convenience, we 
then scale L = \/2Nh so that the critical point is now at h — 1. Choosing cn = 1, 
(66) then reads 

F N (Vmh) = 



B 



N 



-tt 2 (7V+1/2)(V+1)/12/i 2 crC^_ 



f 2Nh) 



2N 2 +N 



Z^-;S V (2N) 



(67) 



The partition function (30) can be analyzed in terms of the same orthogonal 



polynomials appearing in (52). Thus one has [TT] eq. (23)] 



dA 2 



log Z{A; Sp(2A0) 



1 d(+) p(+) 

4/V2 ' N a N+i- 



(68) 



Let us first consider the behavior h > 1 (corresponding to the weak coupling 
phase A < ix 2 ). Knowledge of the large N form of the polynomials for A < n 2 
allows the derivation of the expansion [HI eq. (30)] 



loge -A(V + l/2)(7V+l)/12 z(A . Sp (2JV)) 

A t 1 



-N 



+ 



lir 2 y/irN 



1 



V 1 ~ 



1 



N 



log A 

e -*$hiA/*>) + o (e- 



(69) 



up to terms independent of A (the latter depend on the precise form of cn in (30)), 
where 7(2;) is given by (48). 



us to conclude that for h > 1 



Substituting (69) in (67) we see the first term in the former cancels, allowing 

(70) 



F N (V2Nh) ~ 1 



and furthermore 



'2Nh) ~ e 



-2JV/i 2 7(l/fi 2 ) 



(71) 



where 7(2;) is given in (48). Note that in Model II, the derivative F' N (L) has the 
interpretation of the probability density of the maximal height if at since Fn(L) = 
Pt(Hn < L) is the cumulative distribution of height H^. It is useful now to 
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compare (71) and the analogous result (49) in Model I. Note that, unlike in (49) 



the derivative does not oscillate in sign as N varies. This is reassuring since in the 
present case, F' N (L) has the meaning of a probability density which is necessarily 
positive. 

We turn now to the large deviation formula for h < 1. To derive this large 
deviation formula in this case, in principle we need to repeat, for Sp(2N), the 
strong coupling calculation of Douglas and Kazakov [15J originally done for the 
U(iV) case. In practice, however, as we show below, one can avoid repeating this 
calculation by noting an analogy in the Coloumb gas representation in the two 
cases and thereby relating the behavior of Z{A\ Sp(2iV)) and Z (A; U(N)) in the 
strong coupling regime. As a result, we can then directly use the results of Douglas 



and Kazakov stated in (36) 



Let us first recall that for the ratio of reunion probabilities for Brownian walkers 
on the circle, we relied on the knowledge of the leading asymptotic form of the 
Yang-Mills partition function Z(A; U(iV)) for A > it 2 known from [15]. The latter 
in turn is deduced using the contrained continuum saddle point formula in the 
Coloumb gas representation 

lim ^- 2 \oge A ^ 2 -^Z(A;U(N)) 

= — — J x 2 p(x)dx + J dx J dy p(x)p(y) log \x — y\, (72) 

where the charge density p(x) maximizes the rhs subject to the constraints 

J p(x) dx = l, < p(x) < 1, (73) 

the latter being a direct signature of the spacing in the lattice gas being 1/N (recall 



the discussion above (36)). 



In the case of (30), we begin by supposing all n, > 0, which simply alters 



cat. Then, with rii/(2N) regarded as the continuous variable, the analogue of (72) 
reads 

lim \oge~ A{ - N+l ' 2 ^ N+l ^ l2 Z{A- Sp(2JV)) 

rd rd rd 

= —A I x 2 p(x)dx+ dx dy p(x)p{y) log \x 2 — y 2 \, (74) 
Jo Jo Jo 

where the density p(x) maximizes the rhs subject to the constraints 

d l 
p{x)dx = l, 0<p(x)<-. (75) 
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If we take p(x) = p(—x), p(x) = 2p(x) and d = c then the rhs of (74) can be 
rewritten 



/C PC PC 

x 2 p(x) dx + 2 dx dy p(x)p(y) log \x — y\ 



(76) 



subject to (73). Since the symmetry of the problem implies p(x) is even in (72), 
this is just twice the rhs of (72), so we conclude 

lim log e -MN+i/2)(N + i)/i2 Z{A . Sp(2jV)) 



2 lim -^loge^-^ZiA-MN)). 

TV— >oo iV 



(77) 



Using a similar argument an inter-relation of this type was first noted in |lUj . 
although there it is claimed in Eq. (3.7) that A must be replaced by A/2 on the 
rhs; this in turn is contradicted by the later paper of the same authors |1 1 1 . We 



remark that for A < tt (77) is consistent with (69) and the appropriate case of 

(pl). 



For h < 1, it follows from (77), (36) in the case A> tt 2 , and (67) that 
F N (V2Nh) 



,-2N 2 (F-(l/h 2 )-F+(l/h 2 )) 



(7S 



4.2 Double scaling limit 



Having obtained the large N asymptotic behavior of Fn(L) as a function of L 
in the left (L < y/2N) and the right (L > y/2N) tails, let us now focus on the 
behvaior of F^{L) in the vicinity of the critical point, i.e., when L is close to y2N. 
As in the U(iV) case, this corresponds to the double scaling regime near the critical 



point A 



7T 



in the Sp(2iV) gauge theory. Using results from the gauge theory, 



we will show here in a narrow region \L — \/2N\ ~ N around the critical point 
L = V2N, Fn(L) has the following scaling behavior 



F N (L) -)• T x (^2 11/6 iV 1/6 |L 



(79) 



where the scaling function JFi{t) is precisely the Tracy- Widom distribution function 



for GOE random matrices defined in (12). 



Using results from [15], the double scaling limit of (68) has been analyzed in 
[TT] . In particular, with 



X2N 



n 



2/3 



2N\ 
n c J 



2Ntt 2 
A 



50) 
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one has from Eq. (34) of [TT] that in the double scaling limit N — > oo, x 2 n fixed 



dA 2 



log Z(A- Sp(2iV)) 



71 ■ 



16vr 4 iV 2 



2x 2N k a 



2/3 3/2 

n c 2n c 



fl^N)-^!^))- ( 81 ) 

Tin 



To make use of (81), we first note that with the definition 
W N (A) :=loge-^ N+1 ^ N+1 ^ 12 Z(A; Sp(2iV)) : 



we have from (17) 
Suppose now we fix 



2N 2 + N 12Nir' 
+ 



L 4 



W'J 



2N^ 
V M 2 



16JV 2 7r 4 T „„ /2Nn 2 



L 4 / 3 (l- 



2N\ 



X 2 N ■ 



(82) 



(83) 



4) 



L 2 J ' 

For large N, upon writing the second factor as the difference of two squares and 
setting L = \/2N where this does not lead to a zero term, this is equivalent to 
setting 

X2N = 2 7/6 N 1/6 (L - V2N) . (85) 
'2Nh) 2 in Q, it follows from (ED), Q and 



Recalling (67), and with n c = 
that 

d 2 , fi /rN 4iV 2 + 3iV 
log i^(L = 



<9L 



+ ^ 2 (l - 2^) - L 2/3 (^/f^) + n 2 f[(x 2N )) 



But use of (85) shows that 

AN 2 + W 



+ L 2 (l-2^)=0(N^) 



(86) 



»7) 



allowing these terms to be ignored in ( 86 ) , and leaving us with 
92 log F N (L) = -L 2 ^(^f 2 (x 2N ) + n 2 f[(x 2N ) 



dV 



Finally, use (85) to replace the derivative with respect to L on the lhs of (88) 
by a derivative with respect to x 2 n, and substitute for x 2 n = x according to (57), 
and fi(x) for u(t) = q(t) also according to (57). We then see that (88) implies 

d 2 



dt 2 



logF N ( V2iV(l + t/(2 7 / 3 iV 2 / 3 )) 



1 



q 2 (t) + q'(t) 



■9) 
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Comparison with (11) shows that the distribution of the scaled largest eigenvalue 
in the GOE satisfies the same relation, and so we have 



^iV(l + t/(2 ? / 3 iV 2 / 3 )) ) = F x (t). 



(90) 



As in deducing (60) from (59), to deduce (90) from (89) we have used also the fact 



that the lhs must tend to 1 as t — > oo, which is a consequence of (70). 



We saw in the case of the ratio of return probabilities for the walkers on the 
circle that the large deviation formula relating the values smaller than the mean is 
connected to the left tail of the double scaling distribution about the mean. The 
present problem of the cumulative distribution for the maximum displacement of 
non-intersecting Brownian walkers near a wall exhibits the same feature. Thus 



making use of the expansion ( 38 ) in ( 78 ) shows 



N 



/ 27V(l + t/(2 7 / 3 AT 2 / 3 )) 



e 24 



t— >— oo 



(91) 



which is indeed the leading order tail form of -Fi(t) [50j HI]. On the other hand, 
the right tail of the distribution F^(L) for L > \/2N exhibits a crossover exactly 
similar to the one described above for non-intersecting Brownian motions on a 



circle (63, 65) [55]. Note however that in that case, the second term as in Eq. 
(63) does not oscillate with N because in that case F^(L) has an interpretation 



in terms of cumulative distribution. 



The result (90) is in keeping with known results about fluctuating interfaces 



belonging to the universality class of the Kardar-Parisi-Zhang (KPZ) equation in 
1 + 1 dimensions. Indeed, the top path of such watermelons configuration (see 
Fig. [2]) can be mapped, in the limit N — > 00, onto the height field of such KPZ 
interface in the so-called "droplet" (i.e. curved) geometry [41J. The extreme value 
statistics of such interface in the KPZ universality class and in curved geometry 
has recently attracted some attention [2"B"t EE3"] . In particular, using the fact that 
the maximal value of the height field in the droplet geometry can be mapped onto 
the height field (at a given point) in the flat geometry [33] it was shown, albeit 
indirectly in Ref. [28], that the distribution of (see Fig. [2]), correctly shifted 
and scaled, is indeed described by J-\(t). Here, we obtain this result by a direct 
computation of the distribution of in the large N limit. Moreover the J-\(t) 
fluctuations have previously been established in the case of the distribution of the 
displacement M of the right-most walker amongst N returning vicious walkers in 
discrete time and in the presence of a wall, and with the technical requirement 
that twice the number of walkers be greater than the total number of steps [I]. 
This model gives rise the matrix integral 



(tiS 2N ) 

\ I 5eSp(2A/) 



(92) 
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which upon Poissonization in N is the symplectic analogue of ([TJ. In fact the 
Poissonized form of (92) appears in Hammersley model of directed last passage 
percolation as revised in Section 1.1, but with the points confined to be below the 
line y = x in the square [5], [HI §10.7.1]. 



5 Conclusion 

To conclude, we have studied the normalized reunion probability of N non- inter- 
secting Brownian motions confined on a line segment [0, L\ with three different 
types of boundary conditions : (I) periodic (where the Brownian walkers are thus 



moving on a circle of radius L/2ir) (13), (II) absorbing boundary conditions at 



both extremities and L (16) and (III) reflecting boundary conditions at both 
ends (18). We have shown that, in each of these models, this quantity is given (up 
to a prefactor that we have computed) by the partition function of 2 — d Yang-Mills 
theory on the sphere with a given gauge group, and computed with the heat-kernel 
action. We have found that models I, II and III correspond respectively to the 



group U(iV) ([32]), Sp(2iV) ([33]) and SO(2iV) ([34]). Borrowing results from these 
different field theories, we have shown that, in the large N limit, these reunion 
probabilities exhibit a third-order phase transition as L crosses a critical value 
L c ~ \/N. The region corresponding to L > L c , which corresponds to the weak 
coupling regime in the Yang-Mills theory, describes the right tail of this normalized 
reunion probability, while the region L < L c , which corresponds to the strong cou- 
pling regime, describes its left tail (Fig. [4J. In the critical region of width iV -1 / 6 , 
close to L c , one finds that the reunion probability, correctly shifted and scaled, 
converges to the Tracy- Widom distribution corresponding respectively to GUE, 
J-" 2 (t), in model I and to GOE, J"i(t), in model II and III. One of the main achieve- 
ments of this paper to relate the Painleve equation which describes the singularity 
of the free energy in the double scaling limit of these 2 — d Yang-Mills theories 
with the one defining the Tracy- Widom distributions. In the case of model II, the 
normalized reunion probability has the interpretation of the cumulative distribu- 
tion of the maximal height of the corresponding watermelon configuration (Fig. 
[2]) . Our results thus show directly that this cumulative distribution is given in the 
large N limit by J"i(t), a result was obtained before in a rather indirect way in 
Ref. (2HJ. 

In this paper, we have thus presented the correspondence between boundary 
conditions (in the vicious walkers models) and gauge groups (in Yang-Mills theories 
in two dimensions on a sphere) and discussed its consequences but the deep reason 
behind it deserves certainly further study. An alternative way to explore these 
connections could be to study the relations between vicious walkers models and 
Chern-Simons theory as in Ref. [251 126[ H8] . Yet another point of view could be 
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to adopt the formulation of these vicious walkers problems in terms of Dyson's 
Brownian motion. Indeed, it can be shown that the propagator of this process, 
expressed as a path integral, is precisely given by the partition function of Yang- 
Mills theory on the sphere with the appropriate gauge group, depending on the 
boundary conditions [38]. We hope that this will stimulate further works along 
these directions. 



Acknowledgments 

We thank Alain Comtet for many useful discussions and a careful reading of the 
manuscript. PJF thanks Bernard Jancovici for the opportunity to again visit 
the Laboratoire de Physique Theorique at Orsay, and so make this collaboration 
possible, and also acknowledges the support of the Australian Research Council. 



A Ratio of reunion probabilities for non-inter- 
secting Brownian motions on a circle 



In this appendix, we derive the formula given in Eq. (14) for the ratio Gn(L) 



Riil) / R 1 ^!) introduced in Eq. (13). We consider N non-intersecting Brownian 
motions on a line segment [0, L] with periodic boundary conditions or equivalently 
on a circle of radius L/2tc, starting in the vicinity of the origin at time r = 0. The 
reunion probability -R£(l) denotes the probability that the walkers return to their 
initial configuration at time r = 1. 

Let us begin with the case of iV free Brownian motions, with a diffusion con- 
stant D = 1/2, on a circle of radius L/2ir. The position of the N walkers are thus 
labelled by the angles #i, ■ ■ ■ , On- We denote by Pn(0i, • • • , On] t\pi, • • • , Pn] 0) 
the probability that the positions of the N walkers are 9i, ■ ■ ■ ,9n at time t, 
given that their positions were pi, • • • , Pn at initial time. It is easy to see that 
Pn(@i, ■ ■ ■ ,@N',t\pi, ■ ■ ■ , Pn', 0) satisfies the Fokker-Planck equation 



k=l " 



N 



P N {B X , ■ ■ ■ , N ; t = 0|pi, • • • , p N ; 0) = JJ 5(0 k - p k ) , (93) 



k=i 



which simply comes from the expression of the bi-dimensional Laplacian in terms 
of polar variables (we recall that the radius of the circle is L/2n and the diffusion 
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coefficient is D = 1/2), together with the constraint of periodicity 



Pn(--- A + 2vr, ••• ;t|pi,-- - ,p N ;0) = Pn(- • • A + 27r, • • • ; - • • ,Ptv;0) , 

V 1 < k < N , (94) 

and similarly for a shift of 2n of the variables p&'s. Therefore, from Eq. ([93]) 



and Eq. (94), can be written as the propagator, in imaginary time, of N 
independent quantum free particles on a circle of circumference L. Using path 
integral techniques, one thus writes (using the notation 6 = (9%, ■ ■ ■ ,9^)) 

P N {9 U ■ ■ ■ , 9 N - t\ Pl , ■ ■ ■ , p N ; 0) = (Ole-^lp) , (95) 

with Hl = Yl!i=k ^L,k where tiL,k = — ^ has to be understood as the Hamilto- 
nian of a free particle on a circle of circumference L so that the allowed eigenvalues 

of fiL,k are E nk = 2? ^" fc associated to the eigenvectors <p nk {9) = :^ e " lfe6, > n k £ N. 
Therefore one has in that case 

P N {9t, ■ ■ ■ , N ; t\ Pl , ■ ■ ■ , p N - 0) = *e{9i, • • • , 9 N )V* E (p u • • • , p N ) e - Et , (96) 

E 

with E = E ni + • • • + E nN and ^^(^i, ■ • • , 9^) = (0\E) is the manybody eigen- 
function of Hl- For independent Brownian motions (without the non-crossing 
condition), ^(^i, • • • , 9 N ) is simply the product of the single particle wave func- 

tion • • • ,0*)=n£iW0i)- 

We can now consider the problem of N non-intersecting Brownian motions on 
a circle of circumference L and study the corresponding propagator Pn{0', t\p;0). 



It satisfies the same equations as before (93, 94) together with the additional 
non-crossing constraint: 

P N (9 1 , ■■■ ,9 N ; t\p t , ■ ■ ■ , p N ; 0) = if 9j = 9 k for any pair j, k . (97) 

Following Ref . [47J HQ] , this propagator can be computed using the path- integral 



formalism explained above (96) where, to incorporate the non-colliding condition, 
the many-body wave function ^>e{9i, • • • , On) must be Fermionic, i.e. it vanishes 
if any of the two coordinates are equal. This anti-symmetric wave function is 
thus constructed from the one particle wave functions <p ni of hi by forming the 
associated Slater determinant. Therefore one has, in that case 

1 2tt 2 - 

Me) = 7W.Jt<N^ i6khE = lX^ nl - m 

' k=l 

From the propagator P/v, we can now compute the reunion probability, which 
is the probability that all the iV walkers start and end at the same position on 
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the circle, say Q\ = 6 2 = ■ ■ ■ = On — 0, on the unit time interval. However, such 
a probability is ill defined for a system in continuous space and time. We can go 
around this problem by assuming that the starting and finishing positions (angles) 
of the N walkers are 0<ei<e2<---<eAr and only at the end take the limit 
£j — > 0. Therefore we can compute the ratio of reunion probabilities Gn(L) as 



Gn(L) = lim 



(e|e- Hi |e) 



(e\e- H °°\e) 



(99) 



where denotes the Hamiltonian of the N walkers on the full real axis. Using 



the expressions in Eqs (96, 98), one checks that, in the limit — > 0, powers of 



ti's cancel between the numerator and the denominator in Eq. (99) yielding the 
expression for Gn{L) given in Eq. (14) in the text. 
For instance, for N = 1 one has 



27r r, 2 

e i 2 



(100) 



which can also be written, using the Poisson summation formula 

oo 

G 1 (L) = ^ " 2 • 

n=— oo 

In particular, we obtain the large L behavior as 

G 1 (L) = l + 2e~^ +0(e- 2L2 ) , 



(101) 



(102) 



showing explicitly that G\(L) does not have the meaning of a cumulative distri- 
bution. 



B Reunion probability with absorbing and re- 
flecting boundary conditions 



Here we briefly outline the derivations of the results for Fjy(L) in (17) and En(L) 



in (19). In the first case, we again have N non-intersecting Brownian motions 
on the line segment [0,L] with absorbing boundary conditions at and L. The 
walkers all start at time r = in the vicinity of the origin and we want to com- 
pute their reunion probability -R£ 7 (l) near the origin after time r = 1. The 
calculation proceeds in the same way as in the periodic case in the previous 
appendix. One writes the Fokker-Planck equation for the probability density 
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Pjsr(xi, x 2 , ■ ■ ■ , x^', t\y±, . . . , 2/jv; 0) of reaching {xi, x 2 , . . . , x^} at time t starting 
from the initial positions {y±, y 2 , ■ ■ ■ , Vn}- This reads 



-T N — ~ > , TT^!N 



-T, 

fe=l 

Pat(xi,-- - ,x N ;t = 0\yu-~ ,y N ;0) = Y\_S(x k -y k ). (103) 



k=l 



nk (x) = \l T sin(^) , (104) 



One then writes the solution using path integrals exactly as in (96 ). The rest of the 
calculation is similar as in Appendix A. The only difference is that in constructing 
the Slater determinant, one now has to use the normalized single particle wave 
function as 

2 (UkTlX 

L Sm \TL 

which satisfies the absorbing boundary condition <t> nk (x) = at x = and x = L. 
Using this, one just repeats the calculation of Appendix A to derive the result for 
F N {L) in 0. 

In the reflecting case, the only change is again in the normalized single particle 
wave function that reads 

«U*) = yf cos(^) (105) 

and satisfies the reflecting boundary condition d x <fr nk (x) = at x = and x = L. 
Repeating the rest of the calculation with this modification as in Appendix A, one 
easily derives the result for En(L) in (19). 
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